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SOURCE WAVEFORMS FOR 
ELECTROSEISMIC EXPLORATION 

This application claims the benefit of U.S. Provisional 
Application No. 60/191.041 filed on Mar. 21. 2000. 5 

FIELD OF THE INVENTION 

Tli is invention relates to the field of geophysical pros- 
pecting. More particularly, the invention pertains to source 
waveforms for use in electroseismic exploration. 

BACKGROUND OF THE INVENTION 

The electroseismic method is a geophysical prospecting 
tool aimed at creating images of subsurface formations using 
conversions between electromagnetic and seismic energy. 15 
The electroseismic method is described in U.S. Pat. No. 
5,877,995 (Thompson, et al.). The essence of the electroseis- 
mic method is that high levels of electrical energy are 
transmitted into the ground at or near the surface, and the 
electrical energy is converted to seismic energy by the 20 
interaction of underground fluids, including hydrocarbons, 
with the rock matrix. The seismic waves are detected at or 
near the surface by seismic receivers. To be effective, this 
method requires an input current source with the following 
characteristics: 
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The source should produce large current levels over 
extended time. 



The source should have high electrical efficiency. 

The source should contain little or no DC to avoid plating 30 
the electrode array. 

The frequency content of the source should be adequate 
for the exploration needs. 

The correlation of the source waveform with its reference 
should have sufficiently low side lobe levels. 35 

Little has been published to date on electroseismic wave- 
forms because of the newness of the technique. However, in 
conventional seismic exploration, a seismic vibrator is 
sometimes used as an energy source to generate a controlled 
wavetrain (known as a sweep) which is injected into the 40 
earth. When the resulting recorded seismic data are corre- 
lated with the sweep wavetrain or other reference, the 
correlated record resembles a conventional seismic record 
such as that which results from an impulsive source. 

When a source waveform is correlated with its associated 45 
reference, there will typically be a large peak at the onset 
time of the waveform surrounded by lower peaks at earlier 
and later times. These lower peaks are the correlation side 
lobes. Correlation side lobes are undesirable because they 
can mask smaller desired seismic returns. 50 

It should be noted that the source waveform is just one 
piece of the electroseismic system. Other factors of impor- 
tance include the power waveform synthesizer (that creates 
the waveform) as well as the input electrode array, the 
seismic receiver arrays and various field implementation 55 
issues. 

As stated above, there is litde current technology on 
electroseismic waveforms because of the newness of the 
technique. Some obvious approaches might include pulsing 
or pseudo-random square-wave sequences. Repeated pulses 60 
are inefficient (in energy/second) compared to continuous 
waveforms. Square-wave waveforms would be expensive to 
implement at the required high current levels and would also 
dissipate energy in unwanted high-frequency components. 

What is needed is a source waveform that satisfies the 65 
requirements stated above. The present invention satisfies 
this need. 



In one embodiment, the present invention is a method for 
electroseismic prospecting comprising the steps of (a) 
selecting a source waveform and corresponding reference 
waveform, both chosen to reduce correlation side lobe 
amplitudes, (b) generating an electrical signal based on the 
source waveform, (c) transmitting the electrical signal into 
the ground, (d) delecting and recording the seismic signals 
resulting from conversion of the electrical signal to seismic 
energy in subterranean formations, and (e) correlating the 
resulting seismic signals with a reference waveform to 
produce a correlated seismic record. Preferably, the refer- 
ence waveform is chosen to substantially minimize side 
lobes when correlated with the particular source waveform 
used. 

The source waveform may be constructed from individual 
cycles of a 60 cycles/sec (Hz) sine wave, i.e., standard AC 
electrical power, with the polarity of some such cycles 
inverted as governed by a binary sequence code. The binary 
code is selected to custom design an extended, but finite, 
source wave that has a reference wave that substantially 
minimizes correlation side lobes when the source wave and 
the reference wave are correlated together. The reference 
wave may be the source wave itself or a waveform derived 
from the source wave. Where deeper penetration of the 
subsurface is desired, another embodiment of the present 
invention constructs frequencies lower than 60 Hz by 
switching between the three phases of a 3 -phase power 
source. 

In some embodiments of the invention, the binary 
sequence used is a maximal length shift-register sequence, 
and circular correlation (defined below) is used for the last 
step. In other embodiments, two source waves are transmit- 
ted into the ground. One is a 60 Hz sinusoid wave element 
sequenced by one member of a Golay complementary 
sequence pair; the other is the same wave element sequenced 
by the other Golay pair member. The resulting seismic 
returns are separately correlated with their respective input 
waves and then summed. Both the Golay sequences and the 
maximal length shift-register sequence have excellent cor- 
relation side lobe reduction properties, with the side lobes 
being theoretically reduced to zero in the case of the Golay 
sequences. 

According to the present invention, side lobes may be 
further reduced, for any pseudo-random binary sequence- 
generated source waveform, by making the source wave as 
long as possible before recording equipment limitations 
require that the source generation must be interrupted. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention and its advantages will be better 
understood by referring to the following detailed description 
and the attached drawings in which: 

FIG. 1A illustrates a 60 Hz waveform element, and 
FIG. IB illustrates the autocorrelation of such waveform 
element. 

FIG. 1C shows a waveform element of frequency less 
than 60 Hz, constructed from three-phase 60 Hz waves.- 

FIG. ID shows the autocorrelation of the waveform 
element of FIG. 1C. 

FIG. 2A and FIG. 2B show a Golay complementary 
sequence pair with a 60 Hz waveform element. ^ 

FIG. 3A and FIG. 3B show the autocorrelation of the 
waves in FIG. 2A and FIG. 2B, respectively. 

FIG. 3C illustrates the sum of these autocorrelations. 
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FIG. 4 is a schematic diagram of a shift register. 

FIG. 5A illustrates a maximal length shift-register 
sequence with 60 Hz waveform element, and FIG. 5B shows 
its autocorrelation. 

FIG. 6 A illustrates a modified shift-register sequence with 5 
60 Hz waveform element, and FIG. 6B shows the cross 
correlation of the wave of FIG. 6A with the wave of FIG. 
5A. 

FIG. 7 illustrates a typical field setup for the present 
invention. 10 
FIG. 8 illustrates test results. 

The invention will be described in connection with its 
preferred embodiments. However, to the extent that the 
following detailed description is specific to a particular 
embodiment or a particular use of the invention, this is 15 
intended to be illustrative only, and is not to be construed as 
limiting the scope of the invention. On the contrary, it is 
intended to cover all alternatives, modifications and equiva- 
lents that may be included within the spirit and scope of the 
invention, as defined by the appended claims. 20 

DETAILED DESCRIPTION OF THE 
PREFERRED EMBODIMENTS 

The present invention is a method for determining source 
waveforms to use for electroseismic prospecting for oil and 25 
gas. The source signals taught by the present invention in 
some of its embodiments are in the class of binary-coded 
waveforms. A binary-coded waveform consists of a 
sequence of elements. The individual elements might be, for 
example, a single cycle of a 60 Hz sine wave. In fact, 30 
waveforms that are made up of segments of 60 Hz sinusoids 
(or whatever the local commercial power frequency may be) 
are particularly economical for the electroseismic case 
because this source type can be formed using simple switch- 
ing of commercial powerline signals. These waveform seg- 35 
ments are pieced together with polarities specified by a 
binary sequence. Preferably the binary sequence is designed, 
as explained below, to give minimal side lobes while the 
waveform element is designed to optimize the frequency 
content of the source. 40 

The autocorrelation of a binary-coded waveform will give 
a main wavelet that is the autocorrelation of the individual 
waveform element. (See the discussion below of Golay 
sequences for an example of this.) FIG. 1A shows a single 
cycle 10 of a 60 Hz sinusoid. Its autocorrelation 20 is given 45 
in FIG. IB. This is easily understood by referring to, for 
example, Seismic Data Processing, Ozdogan Yilmaz, Soci- 
ety of Exploration Geophysicists, 1987, 18-19. This wave- 
form element 10 is probably adequate for relatively shallow 
targets. Different elements with lower frequencies can be 50 
constructed when three-phase power is available. An 
example 30 is given in FIG. 1C with its associated autocor- 
relation 40 in FIG. ID. 

By way of further explanation of FIG. 1C, three-phase 
power provides six sinusoids with 60 degrees of phase shift 55 
between them. A variety of approximately square-wave 
signals can be constructed by switching from one sinusoid to 
the next at the crossover point. As one signal starts to fall off 
from its peak, the output is switched to the next signal which 
is rising to its peak. In this manner, an approximate square 60 
wave can be constructed. The square wave can be made with 
a desired width that has an integer number of such cycle 
switches, and hence with a corresponding frequency less 
than 60 Hz. FIG. 1C illustrates an example where the peak 
is prolonged by switching five times to the next-in -phase 65 
sinusoid to yield a square wave with frequency of about 20 
Hz. 



The construction of the waveform element is an important 
aspect of the design of the electroseismic source. Methods 
such as genetic algorithms can be used to determine a 
desirable element for a given target with specified seismic 
attenuation and electromagnetic skin depth. In general, the 
deeper the target, the lower the preferred frequency because 
higher frequencies tend to be absorbed, reducing efficiency. 
Fortuitously. 60 Hz gives good results for targets between 
approximately 100 and 500 feet in depth, for typical sedi- 
ments. Furthermore, the 60 Hz wave element, although not 
preferred, may be used successfully to much greater depths, 
on the order of 5.000 feet: The practical and convenience 
advantages of constructing the waveform from 60 Hz line 
power will be obvious to those trained in the seismic art. 
Furthermore, such hardware implementation is easiest for a 
single frequency sinusoid wave element, in part because 
efficiency does not have to be sacrificed for broadband 
amplification. 

Correlation side lobes are of critical importance in elec- 
troseismic exploration because there can be a very large 
peak at zero time. This large peak comes from unavoidable 
direct pickup at the receivers from fields related to the input 
currents. The large peaks will have significant correlation 
side lobes. Even though these direct-pickup side lobes are 
reduced from the peak amplitude, they may still be large 
enough to mask the much smaller desired electroseismic 
returns. The level of the direct pickup can be moderated by 
proper field design and/or by other innovations such as 
receiver modifications; nonetheless, it is best to minimize 
the impact of the direct pickup by using a source waveform 
with minimal correlation side lobes. (Correlating with the 
appropriate reference waveform for the source waveform 
will reduce all side lobes, including the direct pickup side 
lobes because the direct pickup also is caused by the applied 
signal.) 

In the preceding discussion of correlating a source wave 
with a reference wave to minimize side lobes, the source 
wave is a surrogate for the electroseismic response. In the 
present invention, it is the recorded seismic response to the 
electrical source wave that is correlated, in a later processing 
step, with the chosen reference wave. The present invention 
is based on detecting the linear component of the electroseis- 
mic response. As a consequence of linearity, such response 
will be proportional to the source wave. Therefore, custom- 
designing a source wave and a reference wave to have a 
large central peak and minimal side lobes when correlated 
together insures that the seismic response to that same 
source wave will similarly produce a large central peak and 
minimal side lobes when correlated with that same reference 
wave. 

There are at least two types of binary sequences that 
according to the present invention are well suited for elec- 
troseismic waveforms. One is the maximal length shift- 
register sequence which is found in Gdlomb, S., Digital 
Communications with Space Applications, Prentice Hall, 
Inc. (1964). The other is the Golay complementary sequence 
pair, which is found in Golay, M. J. E., "Complementary 
Series", IRE Transactions on Information Theory (1961) 
Vol. 7, 82-87. These sequences can be configured to give 
substantially minimal correlation side lobes. 
Golay Sequence 

In one embodiment, the present invention uses sinusoidal 
wave elements sequenced by Golay complementary 
sequence pairs. A complimentary pair of Golay series is 
defined by Golay as a pair of equally long, finite sequences 
of two kinds of elements which have the property that the 
number of pairs of like elements with any one given sepa- 
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ration in one series (viewing each series as cyclical for 
purposes of determining separation) is equal to the number 
of pairs of unlike elements with the same given separation 
in the other series. These sequence pairs have the property 
that the sums of the sequence autocorrelations have zero side 5 
lobes. One such pair of length 8 is {-1 -1 -1 1-1-1 1-1} 
and {-1-1-1 1 1 1-1 1}. This pair of series can be seen 
to satisfy the definition given above. For example, taking a 
separation of 3 spaces, the first series has 6 like pairs and 2 
unlike pairs while the second series has 2 like pairs and 6 10 
unlike, and so on. This pair of series of length 8 is shown in 
FIGS. 2A and 2B with a 60 Hz element. FIG. 2A shows the 
first-mentioned sequence 50; FIG. 2B shows the second 
mentioned sequence 60. A "1" in either sequence means a 
single cycle of 60 Hz sinusoid with normal polarity, illus- 15 
trated by the cycle starting at 52 and ending at 54 in FIG. 2A; 
a "-1" means a 60 Hz wave element with inverted polarity, 
illustrated by the cycle starting at 54 and ending at 56 in FIG. 
2A. 

The autocorrelations and sum are shown in FIGS. 3A-C. 20 
FIG. 3 A shows the autocorrelation 70 of the extended 
waveform segment 50; autocorrelation 80 of extended wave- 
form segment 60 is shown in FIG. 3B. FIG. 3C depicts the 
sum 90 of autocorrelations 70 and 80. In theory, the side 
lobes cancel completely. The sequence-pair approach is also 25 
very efficient in the sense that the waveform is continuous, 
i.e., without breaks as contrasted with, for example, the 
modified shift-register sequence described below. The draw- 
back is that the side lobe cancellation is theoretical, depend- 
ing as it does on proper subtraction of the relatively large 30 
side lobes of autocorrelations 70 and 80. This can be 
problematic in practice where there can be frequency drifts 
in the power supply or fluctuations in the signal amplitude. 

Registering the sequence to a reference line frequency, 
scaling the data to match amplitudes, and selecting particu- 35 
lar Golay sequence pairs that have relatively low side lobes 
at the lags of interest can minimize these problems. 

For a typical field implementation, a preferred embodi- 
ment of the present invention might include a 60 Hz wave- 
form element and a Golay sequence of length 1664 (27.7 40 
seconds). Golay sequences can be constructed of smaller 
sequences using methods discussed in Kounias, S., 
Koukouvinos, C., and Sotirakoglou, K., "On Golay 
Sequences", Discrete Mathematics, (1991) Elsevier Science 
Publishers B.V., Vol. 92, 177-185. Golay sequences exist 45 
only in certain lengths, e.g., there are Golay sequences of 
length 8 but none of length 6. Golay showed that sequences 
exist for lengths given by 2 / 10* 26\ where j, k, and 1 are 
nonnegative integers, i.e., 0, 1,2... (see Kounias, p. 178). 
Other lengths that can, but not necessarily do, have Golay 50 
sequence pairs are given by a 2 +b 2 where a and b are integers 
(including zero). For lengths up to 50, Golay claims 
sequence pairs exist for the following lengths: 2, 4, 8, 10, 16, 
1 8, 20, 26, 32, 34, 36, 40 and 50 (see Golay, pg. 84). Kounias 
provides an algorithm for generating solutions for a given 55 
length (see Kounias, pg. 184). Lemma 1 in the Kounias 
paper (pg. 178) shows how to construct longer Golay 
sequences from shorter sequences. A Golay sequence pair of 
length nm can be constructed from a Golay pair of length n 
and another Golay pair of length m. Following these con- 60 
sanction rules, there are many possible sequence pairs of a 
specified length. 

A particular pair for a given length can then be selected 
using, for example, an exhaustive search to select a pair with 
minimal side lobes prior to cancellation. This selection will 65 
nunimize the residual side lobe energy from imperfect 
cancellation. Each of the Golay-pair sequences is run in the 



field as a separate source with some dead time allotted after 
each to collect seismic returns. The reference correlations 
and summing are processing steps. 

As an example of how to generate other complementary 
Golay series of a given length, consider the pair of length 8 
given above: 

-1 -1 -1 1-1-11 -I and -l -i -l | | | -i i 
Golay showed that any of the following operations on a 
given pair produces another pair of complementary series: 

(a) Interchanging the series: 

(b) Reversing the first series: 

(c) Reversing the second series: 

(d) Altering (replacing each element by its opposite) the 
first series: 

(e) Altering the second series: or 

(f) Altering the elements of even order in each series. 
Following these rules, six more complementary series of 
length 8 are: 

-1-1-1111-11 and -1 -1 -I 1 -1-1 1 -1 
-1 1-1-11 -1-1 -land -1-1-1111-11 
-1 -1 -1 1 -1 -1 1 -1 and 1 -1 1 1 1 -1 -1 -1 
1 1 I -1 1 1 -1 1 and -1 -1-1 111-11 
-1 -1 -1 1-1 -1 1 -1 and 1 1 1 ^1 -1 -1 1 -1 
-11-1-1-1111 and -111-1-11 -1 -1 -1 
Many more complementary series of length 8 (some of 
which may be identical ) can be generated by performing two 
or more of the six operations given above. In addition, these 
operations can be performed on the sub-sequences (e.g., 
pairs of length 4 and 2) prior to combination to make a 
length 8 sequence. 

Maximal Length Shift-Register Sequence 

Golomb defines a pseudo-random binary sequence 
("PRBS") as any binary sequence generated by a determin- 
istic process (such as a shift register) in such a way that the 
sequence will satisfy whatever tests for randomness that 
may be selected (see Golomb, pp. 7-16). A shift register of 
degree n is a device consisting of n consecutive binary (1,-1 
or 1, 0) storage positions or "registers", which shifts the 
contents of each register to the next register down the line, 
in time to the regular beat of a clock or other timing device. 
In order to prevent the shift register from emptying by the 
end of n clock pulses, a "feedback term" may be compiled 
as a logical (i.e., Boolean) function of the contents of the n 
positions and fed back into the first position of the shift 
register. 

For example, consider the case where n=4 and the feed- 
back function is to add the contents of the third and fourth 
registers, the sum to become what is put into register 1 after 
the next shift empties it. Such addition of binary numbers is 
called modulo 2 addition and is denoted by the symbol 0. 

Thus in the binary {1,0} domain, 000=0; 001=100=1; 
and 101=0. Such a shift register is illustrated in FIG. 4. 

It can be shown that this feedback function can be 
expressed as the following recursion formula: 

where X, is the contents of any one of the four registers for 
the i-th shift. Thus, the contents of any register are the 
modulo 2 sum of what was in that same register three shifts 
previously and what was in that same register four shifts 
previously. 

Starting the process with the contents of all four registers 
set to 1, i.e., X c (X l )=X 0 (R 2 )=X o (R3)=X o (X 4 )=l, the four 
registers take on the following values before the numbers 
begin repeating: 
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The numbers generated in register 1 (the other registers 
generate the same sequence with cyclic permutation) are the 
"shift-register sequence" for this particular shift register of 
length 4 and particular recursion relationship. It can be seen 
that the number in register 1 for any value of i is the modulo 
2 sum of the numbers in register 3 and 4 one shift previously 
which, in turn, are the same two numbers that were in 
register 1 three and four shifts previously, as required by the 
recursion formula. This sequence satisfies, for example, the 
following three randomness tests which are proposed by 
Golomb at page 10 of his book: 

1. In each period of the sequence, the number of ONE's 
differs from the number of ZERO'S by at most 1. 

2. Among the "runs" of ONE's and of ZERO'S in each 
period, one-half the runs of each kind are of length one, 
one-fourth are of length two, one-eighth are of length 
three, etc., as long as these fractions give meaningful 
numbers of runs. 

3. If a period of the sequence is compared, term by term, 
with any cyclic shift of itself, the number of agreements 
differs from the number of disagreements by at most 1. 

In the sequence of length 15 generated above (register 1), 
there are 8 ONE's and 7 ZERO'S, satisfying test 1. The 
sequence has 5 runs of length one, 2 of length two, and 2 of 
length three, which closely satisfies test 2. (A "run" occurs 
where a number is repeated in the sequence, except for a run 
of length one which is a single occurrence.) Comparing the 
sequence of register 1 with that of register 4, which is a 
cyclic permutation of the register 1 sequence, shifted 3 
positions, one finds 7 agreements and 8 disagreements, 
satisfying test 3. This consistent satisfaction of randomness 
tests is a difference between pseudo-random series and truly 
random series. Because they are deterministic, each and 
every pseudo-random series will satisfy the test Truly 
random series will satisfy the tests only on average. 

The output of any shift register is ultimately periodic, with 
a period not exceeding 2 N where n is the degree, or length, 
of the shift register (Golomb, p. 9). For linear recursion 
formulas, defined by Golomb at page 9, the period is at most 
2"-l. In the example above, where n=4, the period is 15 and 
therefore the sequence generated above has the maximum 
possible length, and accordingly is called a maximal length 
shift-register sequence. 

An example of a maximal length shift-register sequence 
oflength7is{-l 1-1 1 1 1 -1}. FIG. 5 A shows the resulting 
extended waveform segment 100 using a 60 Hz element. 
The circular autocorrelation 110 of waveform 100 is shown 
in FIG. 5B. The central portion of waveform 110 is the 
autocorrelation of a 60 Hz cycle and the side lobes are 60 Hz 
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with relative amplitude of 1/7 (for a length 7 sequence). This 
level of side lobes might be acceptable for long sequences 
but alternate approaches are possible. 

Foster and Sloan, for example, altered waveform 100 to 
5 include only the positive binary elements 120. with the 
negative elements replaced by zero-amplitude elements, as 
shown in FIG. 6A. (Foster, M. R.. and Sloan. R. W.. "The 
Use of Pseudonoise Sequences to Code a Pulsed Neutron 
Logging Source". Geophysics (1972) Vol. 37. 481-487). 
When circularly correlated using waveform 100 (FIG. 5A) 
as the reference, the result 130 now has zero side lobes (FIG. 
6B). The disadvantage of this approach is that the peak value 
is roughly halved because of the zero-amplitude elements. 

To explain terms used above, autocorrelation means the 
correlation of a signal with itself. Circular correlation can be 

15 explained as follows: In a standard correlation process, the 
signals are assumed to be "zero-padded" prior to correlation, 
i.e., the signal wave sequence and its reference are assumed 
to drop to zero amplitude before and after the sequence. The 
correlation process involves the cross product of one signal 

20 and a shifted version of the second signal for various shifts. 
With the zero padding, the portion of the shifted signal that 
passes the end of the stationary signal has no effect because 
it is multiplied by the appended zeros. In the case of circular 
correlation, the signals are assumed to repeat rather than 

25 have zero padding. Thus, as the shifted signal passes the end 
of the stationary signal in the correlation process, it begins 
to overlap the beginning of the stationary signal. Where 
circular correlation is used, it is used because it reduces the 
side lobes better than standard correlation. 

30 The choice of a binary sequence and waveform element 
will depend on the target, the geology, and the field con- 
figuration. The decision on a waveform element would be 
based largely on the target depth and the expected attenua- 
tion. The decision on the binary sequence type would be 

35 based on the relative level of direct electromagnetic pickup 
and, consequently, the need for side lobe reduction. An 
appropriate selection will likely require computer modeling 
and field tests. As discussed further below, the longer the 
sequences, the more the side lobes will be reduced. Capa- 

40 bilities of the recording equipment are a practical limit on 
sequence length. Although the present invention is primarily 
for use on land, it can give useful results in a marine" 
environment. 

In general, longer binary sequences tend to produce lower 

45 correlation side lobes than shorter sequences. In the case of 
waveforms generated by maximal length shift-register 
sequences, it is known that the amplitudes of the side lobes 
after the sequence is circularly correlated with itself varies 
as 1/L where L is the length of the sequence. Therefore, 

50 longer maximal length shift-register sequences are pre- 
ferred. Truly random sequences are known to exhibit a lesser 
reduction with length, varying inversely as the square root of 
the length, but this effect is manifested on a statistical basis. 
The other preferred type of binary coding for the present 

55 invention, Golay sequence pairs, are neither pseudo-random 
nor random. They also tend to show an inverse dependence 
of side lobe amplitude with sequence-length, but the depen- 
dence is much weaker than 1/L. Although side lobe cancel- 
lation for Golay pairs is theoretically perfect, regardless of 

60 sequence-length, it is unlikely that this perfect cancellation 
will be achieved in practice as pointed out above. Therefore, 
Golay sequences are preferably chosen to have minimal side 
lobes before the pairs are added together, and sequence 
length is one factor in such a selection. Other than length, the 

65 basis for choosing a sequence is the same as with PRBS 
sequences: one must try various sequences and observe what 
side lobe attenuation is achieved. 



US 6,477,113 B2 



10 



The correlation for Golay sequences is standard correla- 
tion. Pseudo-random sequences need circular correlation to 
best reduce side lobes, which tends to place a practical limit 
on the length of such sequences. A pseudo-random sequence 
needs to be repeated (in the field) at least once to provide 5 
circular overlap. The first cycle can be used only to provide 
delayed correlation effects on the next cycle, and hence is 
lost for data-gathering purposes. While a longer PRBS is 
preferred for side lobe reduction, a shorter PRBS is preferred 
to avoid data loss and waste of field time. The preferred 10 
compromise is to use an intermediate length PRBS and 
repeat it several times, perhaps three to seven cycles in all. 
Thus a maximal-length shift register sequence of length 255 
might be selected and repeated six times for a total of seven 
cycles. This would consume about 28 seconds in the field 15 
assuming a 60 Hz wave element, and only 1/7 of the data 
would be. in effect, lost. The data thus gathered from such 
a PRBS example would compare fairly closely in quantity to 
that obtained from a Golay-pair sequence of length 1664 (as 
6x255=1530 compares to 1664). 20 

FIG. 7 is a diagram of a possible field layout for deploy- 
ment of the present invention. Electric current is injected 
into the subsurface 140 by applying a voltage from a power 
source 141 between two buried electrode wires 142 and 143. 
The electrode wires are typically bare 4/0 copper cable. 25 
Other sizes of cable may be used as needed to conduct the 
needed current. The electrode wires should be buried just 
deep enough to make good electrical contact with the soil. 
Typically, this depth will be in the range from 1 inch to 10 
feet, but in some cases, such as particularly dry surface soil, 30 
even deeper burial may be desirable. 

FIG. 7 shows the current paths 144. The direction shown 
is that of the flow of electrons, from the negative electrode 
to the positive electrode. The power source 141 provides the 
current and consists primarily of a waveform generator 35 
capable of generating the binary-coded waveform segments 
preferred by the present invention. The ultimate power 
source is typically the local electrical utility power lines 
(connection not shown in FIG. 7). One or more generators 
may also be used. The current paths 144 shown in FIG. 7 40 
represent those current paths that penetrate down to the 
depth of the target 145, typically a petroleum reservoir. Such 
current trajectories will be nearly vertical below the elec- 
trode. Experience shows that the maximum vertical current 
146 will typically be directly under or even slightly to the 45 
outside of the electrode wires. Accordingly, the maximum 
amplitude of the resulting surface-directed seismic wave 147 
(similar wave below positive electrode not shown) will 
occur along the line of maximum vertical current, and this 
determines the preferred location of the seismic detectors 50 
148. 

The seismic detectors may be placed anywhere on the 
surface, but the preferred location is outside the electrode 
wires rather than between the electrode wires. The seismic 
detectors may be geophones, hydrophones, accelerometers 55 
or any similar device. Such seismic equipment is well 
known to those skilled in the art. Preferably, the seismic 
detectors are buried beneath the surface to reduce seismic 
noise. 

Typically, the configuration shown in FIG. 7 is designed 60 
to cover the entire area of interest i.e., source and receivers 
will not need to be repeatedly moved to progressively cover 
the area of interest as with conventional seismic. This is one 
advantage of the electroseismic method, although a single 
setup is not essential to using the present invention. (The 65 
electrode wires and receivers may be dug up when the 
experiment is over and moved to other locations.) 



Accordingly, the length of the electrodes may vary between 
one-tenth of the reservoir (target) depth to several times the 
reservoir depth. The separation between electrodes in pre- 
ferred embodiments of the present invention will be approxi- 
mately equal to the target depth. 

In electroseismic surveying, the electromagnetic source 
wave reaches all of the target at essentially the same instant 
of time. One is therefore usually interested only in the 
upward-traveling seismic waves that arrive at all of the 
geophones at approximately the same time, assuming the 
geophones are deployed as in FIG. 7. Thus, in processing, 
the ordinary seismic noise can be filtered out because it 
exhibits what in the seismic art is called "moveout". The 
source-to-receiver distance, called "offset", is small in elec- 
troseismic prospecting. In a conventional seismic survey, a 
much larger surface area would have to be covered because 
the longer offsets would be needed from every shot point to 
provide sufficient signal-to-noise ratio in common mid-point 
gathers. The area reduction in electroseismic vs. conven- 
tional seismic may be approximately fourfold. 

The binary-coded waveforms discussed above meet the 
five requirements for electroseismic exploration stated pre- 
viously. The need for large current levels is addressed since 
these are continuous waveforms instead of, for example, 
pulses that would have significant dead time. The use of 
simple 60 Hz (or constructed 3-phase) elements also allow 
large current levels and high electrical efficiency since the 
related hardware is simplified. The lack of DC is assured - 
since each waveform element (e.g., a full cycle at 60 Hz) has 
no DC component (i.e., its mean value is zero); it follows, 
therefore,, that a set of such elements would have no DC 
component. The frequency content of the source can be 
matched to the exploration target by adjusting the frequency 
of the waveform element. Finally, minimization of side 
lobes has been discussed at length. 

To further explain the preceding statement about adjusting 
the frequency of the source to achieve the desired depth 
penetration, neither the source wave nor the return seismic 
response is composed of a single frequency. The phase 
inversions and (in some embodiments) the zeroing of certain 
elements produce waves composed of many frequencies in 
the sense of their Fourier analysis decomposition. This is 
necessary to the invention, i.e., that the waves have a 
bandwidth of frequencies. If the seismic return wave were a 
single frequency, there would be no wave that would cross- 
correlate with it to produce a localized pulse in the process- 
ing step. A bandwidth of frequencies is needed to produce 
the desired spike. Elementary Fourier analysis teaches that 
the sharper the spike, the wider the needed bandwidth. Thus, 
although the source and return waves each have a spread of 
frequencies, it is reasonable to expect, and Fourier decom- 
position can prove, that the frequency distribution of both 
waves will peak at the frequency of the building block, the 
waveform element. Thus, the desired subsurface penetration 
may be achieved by varying the frequency of the waveform 
element. 

FIG. 8 illustrates the results obtained by applying the 
present inventive method at the Friendswood gas field in 
Texas. The region depicted is in the vicinity of Well No. 181. 
The site was selected based on knowledge of gas leakage at 
the surface indicating shallow gas deposits located above the 
producing formations. A conventional seismic survey of the 
area showed seismic bright reflections shown in FIG. 8, of 
which band 160 is by far the most prominent. The five dark 
bodies 170 represent regions of high amplitude electroseis- 
mic signal. Electroseismic signals of less than about half the 
peak value are made transparent by the display. The detected 



